Anatomical traits explain drought response of seedlings from wet tropical forests

Abstract Water availability regulates plant community dynamics but the drought response of seedlings remains poorly known, despite their vulnerability, especially for the Asian tropics. In particular, discerning how functional traits of seedlings mediate drought response can aid generalizable predictions of tree responses to global environmental change. We assessed interspecific variation in drought response explained by above‐ and below‐ground seedling traits. We conducted a dry‐down experiment in the greenhouse using 16 tree species from the humid forests of Western Ghats in southern India, chosen to represent differences in affinity to conditions of high and low seasonal drought (seasonality affiliation). We compared survival, growth, and photosynthetic performance under drought and well‐watered conditions and assessed the extent to which species' responses were explained by seasonality affiliation and 12 traits of root, stem and leaf. We found that the species from seasonally dry forest reduced photosynthetic rate in drought compared with well‐watered conditions, but seasonality affiliation did not explain differences in growth and survival. Performance in drought vs well‐watered conditions were best explained by anatomical traits of xylem, veins and stomata. Species with larger xylem reduced their growth and photosynthesis to tolerate desiccation. In drought, species with smaller stomata correlated with lower survival even though photosynthetic activity decreased by a larger extent with larger stomata. Overall, anatomical traits of xylem and stomata, directly related to water transport and gas‐exchange, played a more prominent role than commonly used traits (e.g., specific leaf area, leaf dry matter content) in explaining species response to drought, and may offer a good proxy for physiological traits related to drought tolerance of seedlings.


| INTRODUC TI ON
Water availability drives the distribution and performance of plant species, even within relatively wet biomes (Comita & Engelbrecht, 2009;Engelbrecht et al., 2007;Esquivel-Muelbert et al., 2019;Gopal et al., 2023;Krishnadas et al., 2021).Quantifying plant response to water stress is imperative given that global environmental change will alter patterns of precipitation and seasonal drought (Hajek & Knapp, 2022;Mishra et al., 2020).For closed canopy forests, much attention has focused on the response of mature trees to drought (Anderegg et al., 2015;Hartmann et al., 2018;Rowland et al., 2015;Suresh et al., 2010), but less is known about seedlings.Seedlings can be highly vulnerable to drought because their shallow roots do not allow access to deeper water layers during dry periods.Moreover, the light-limited forest understory may impose carbon starvation during drought (Amissah et al., 2015;Slot & Poorter, 2007) such that even short dry spells can result in substantial seedling mortality (Engelbrecht et al., 2007).Seedling responses to water availability can therefore affect regeneration and persistence of species in a changing climate and shape community structure of later life-stages (Comita & Engelbrecht, 2009, 2014).
Interspecific variation in drought response can be mediated by species' traits-anatomical, morphological and physiological characteristics associated with differences in strategies of resource use (Poorter & Markesteijn, 2008;Slot & Poorter, 2007).For example, wood density and leaf mass per area (LMA) are correlated positively with drought tolerance (Muscarella et al., 2016;O'Brien et al., 2017), but may not play a direct function in the plant's water use efficiency or transport, whereas variation in anatomical features of stomata and xylem are functionally related to water transport and use of water and help regulate water balance and response to dry conditions.
Stomatal and xylem traits may thus serve as good proxies for hydraulic function (Bartlett et al., 2016) but are not well-characterized for tropical tree seedlings, especially in terms of seedling response to drought.In addition, the importance of roots in accessing ground water is commonly understood, but the role of root traits in mediating species' drought responses remains poorly explored.Root traits can be a key component of plant functional strategies and wholeplant response to drought stress (Freschet et al., 2021;Laughlin et al., 2021;Weemstra et al., 2023), which may be logistically difficult to study for mature trees but convenient for seedlings.
Evergreen species may have evolved to withstand seasonal droughts either by avoiding or tolerating desiccation, with functionally distinct strategies arising from different trait combinations.
Species from less seasonal habitats may try to maintain water status through a drought to avoid desiccation.These species may maintain photosynthesis and growth and hence adapted to have smaller but more numerous stomata, highly conductive larger xylem to transport adequate water, and deep roots to maintain access to water (Chitra-Tarak et al., 2021;Reich, 2014;Reich et al., 2003).Species from more seasonal forests may have evolved to tolerate desiccation by decreasing photosynthesis to reduce water loss, with smaller stomata, and smaller xylem with lower conductivity (Reich et al., 2003).
Desiccation tolerant evergreen species from more seasonal forest should be resource-conservative to prevent carbon starvation.Which should be reflected in features such as lower investment in below-ground tissue (root mass fraction, RMF) to forage for water compared with desiccation avoiders (Lawson & Blatt, 2014;Rodríguez-Ramírez et al., 2022), as well as other associated traits such as lower specific leaf area (SLA) and specific root length (length per unit mass, SRL), higher leaf dry matter content (LDMC) and stem tissue density (Poorter & Markesteijn, 2008).Alternatively, traits may not be coordinated and different trait combinations may offer equivalent responses to drought, with individual traits varying in their relationship to drought response (Laughlin et al., 2020).
Easily measurable traits relevant to resource-use may therefore help categorize species as avoiding or tolerating desiccation based on a whole-plant evaluation of photosynthesis, growth and survival (Bartlett et al., 2012;Sun et al., 2020).
We assessed the relationship between traits and species' performance in dry vs. well-watered conditions for seedlings of 16 tree species from a wet tropical forest region in a global biodiversity hotspot in southern India, the Western Ghats.Across the region, most seedlings establish during the monsoon, but then have to withstand dry periods varying from 3 to 9 months to persist.Given the range of dry months in this landscape, the seedlings must have a suite of traits evolved to tolerate or avoid desiccation.The species were chosen to represent a gradient in their adaptation to seasonal drought, climatic water deficit during the dry season (henceforth seasonality), based on a previous analysis of tree species distributions in the Western Ghats (Krishnadas et al., 2021).We asked: 1. Do species from less vs. more seasonal forests vary in their drought response?
We expected that species associated with greater seasonality would tolerate desiccation better than species from less seasonal forests, and therefore survive better but decrease photosynthesis and/or growth in drought vs. well-watered conditions.

How well do individual traits explain interspecific variation in drought response?
We expected that species with higher SLA, lower LDMC, higher root mass fraction (RMF) or deeper roots, larger stomata and wider xylem, suite of traits found in plants adapted to wetter habitats, would maintain better function (photosynthesis and growth) early in drought, but survive less than species with lower SLA, higher LDMC, and lower RMF or shallow roots, smaller stomata, and narrower vessels.
3. Do distinct trait combinations characterize species that avoid or tolerate desiccation?
We expect that traits associated with conservative water-use or safe hydraulic strategy, e.g., thicker leaves, smaller stomata and xylem, would correspond with responses characteristic of desiccation tolerant species adapted to more seasonal habitats-decreased photosynthesis and growth.Species adapted to less seasonal habitats, with thinner leaves and larger stomata or xylem would avoid desiccation by functioning at similar levels in drought compared with watered conditions.We also expect that species with higher SLA, lower LDMC, deeper roots and higher RMF would avoid desiccation.Species with lower SLA, higher LDMC, shallow roots and lower RMF were expected to tolerate desiccation.

| ME THODS
We conducted a dry-down experiment in the greenhouse to gauge the response of seedlings to seasonal drought.Hence, our experiment represented a scenario where seedlings growing in a limited volume of soil were exposed to sudden drought and results should be interpreted accordingly.Seeds of 16 different tree species (Table A1) were obtained from two locations, Valparai in southern Western Ghats and Sirsi in northern Western Ghats, representing less and more seasonal climates respectively.Sirsi has higher mean annual rainfall but more negative climatic water deficit (CWD) (Sirsi: Mean annual precipitation (MAP) = 3425 mm/year, number of dry months: 6, CWD = −826.9mm/year; Valparai: MAP = 2169.3mm/ year, number of dry months: 4, CWD = −396.51mm/year, Source: Chave et al., 2014;Fick & Hijmans, 2017), which is a measure of both length and magnitude of drought stress on plants (Stephenson, 1998) and characterizes seasonal drought conditions well (Krishnadas et al., 2021).

| Greenhouse experiment
Between July and November 2019, seeds were germinated in germination trays containing a 2:1 mixture of all-purpose soil and cocopeat at the greenhouse at the National Centre of Biological Sciences (Bengaluru, India).Because species fruited at different times, had recalcitrant seeds, and germination times varied among species, we planted species for germination whenever seeds became available.
Seeds were watered every alternate day and monitored for germination.Once germinants lost their cotyledon leaves and developed at least four mature leaves they were transplanted into PVC pipes of diameter 15 cm and depth of 45 cm (Pot volume = 7.95 L).The PVC columns were filled with a prepared soil mixture (all-purpose soil: sand: cocopeat: compost = 1:1:1:0.1)and randomly assigned to six blocks within the greenhouse.Seedlings within a block were randomly assigned to drought and control treatments.Each treatment per block had an equal number of individuals of a species, but numbers differed among species depending on seedling availability.To maintain a consistent post-germination age at which seedlings were subject to drought, species were included in the experiment as their seedlings reached four-leaf stage.All new seedlings were allowed to acclimate in their pots for 2 weeks and watered every 2 days to avoid transplant shock.Seedlings that perished during this period were replaced and included after acclimation.In total, we had 434 seedlings in the experiment with 1-3 individuals per species per treatment per block.
After the acclimation period, we stopped watering the seedlings in the drought treatment for each block.Seedlings in the control treatment were given 75-100 mL of water twice weekly to ensure a consistent level of soil moisture.The greenhouse was maintained at approximately 65% relative humidity and 26°C temperature.
Initially, we planned the drought treatment to last 16 weeks to simulate the post-monsoon seasonal drought that occurs from January through April in the Western Ghats.However, monitoring and measurements were disrupted due to the COVID-19 lockdown.For 12 of the 16 species, we implemented the experiment for 16 weeks.
For the remaining four species, control plants were watered by the greenhouse staff but data collection was not possible due to restricted access to the facility.We harvested these plants in the 12th week of their drought experiment when we got permission to access the facility.These 12 weeks still give us a reasonable estimate of interspecific variation in response to seasonal drought conditions.
For growth, stem height, leaf length and width (for the four largest leaves) and survival were recorded weekly for the first 8 weeks and then once a month till the plants were harvested.We calculated the growth rate at the end of 12 weeks for all individuals as follows: (final height -initial height)/time in weeks.Seedling final height was strongly correlated with initial height (r = 0.7, t = 20.8,p < .001),hence we calculated relative growth rates, standardized by initial height as follows: (final height -initial height)/(time in weeks*initial height).Photosynthesis was measured monthly using a LiCOR 6400XT on seedlings of all species from control and drought treatments.Photosynthesis rates were measured with a 2 × 3 cm 2 leaf chamber.The chamber conditions were maintained at PAR 1500 μmol m −2 s −1 , 65% humidity, 400 μmol s −1 flow rate, fan speed at 10,000 rpm.Measurements were taken between 0700 and 1200 h.

| Trait measurements
For leaf trait measurements (listed in Table A2), 1-3 leaves each were collected from selected seedlings and immediately stored in a ziplock bag between moist tissue paper to make sure the leaves were fully saturated with water (Cornelissen et al., 2003).Five seedlings from the control group were selected for each species for trait measurements.
To measure SLA (cm 2 g −1 ) and LDMC (μm mg −1 ) (Cornelissen et al., 2003), the fresh leaf area was measured using the CI-202 leaf scanner (CID bio-science).After measuring fresh weight, leaves were placed in paper bags and air-dried for a week, following which they were dried at 60°C for 3 days to get dry leaf mass.SLA was calculated as the ratio of leaf area to its dry biomass.LDMC was calculated as the ratio of dry leaf biomass and fresh leaf biomass.
To measure stomatal traits, the leaf surface was wiped thoroughly to get rid of dust and excess moisture.Nail polish imprints of the underside were prepared and scanned under 40× magnification in a light microscope.For stomatal density (SD, number per cm 2 ), four fields of vision were imaged from each imprint.The number of stomatal apertures was counted and the image area was measured.We chose three of the clearest stomates from each of the images and estimated stomatal size (SS, μm 2 ) by drawing an ellipse around the stomata and stomatal length (SL, μm) was measured along the long axis of each stomata.We then calculated stomatal area fraction (SAF) for cm 2 of leaf by taking a product of SD and SS (SD * SS * 10 −8 ).Model results for SAF are similar to SD (Figure A1) and therefore we decided to use SD for further analysis and interpretation.
To measure vein density (VD, mm −2 ), whole leaves were immersed in petri dishes with cleaning bleach and monitored daily.
Once the leaves lost all the chlorophyll, they were washed with distilled water and stained with methylene blue diluted 1:100 with distilled water.Whole leaf samples were fit on to clean slides and imaged under 20× magnification.Veins of all the orders were measured and density was calculated per mm 2 .Xylem diameter (XD, μm) was measured using an inch-long section of the stem that was cut off and sectioned immediately into thin slices with a stainless-steel blade.
Five to ten of the thinnest sections were observed under the 40× light microscope and the diameter of 8-10 of the largest xylem vessels were measured using the line and circle tools of the software.All microscopic measurements were performed using a Primostar, Carl Zeiss microscope and images were taken using the camera accessory Axiocam-105 mounted on the microscope.All images were analyzed to measure the stomatal, vein and xylem traits using the ZEISS ZEN2 (v1.0, RRID: SCR_013672).
To measure root traits, harvested roots were stored in a freezer until scanning.The root samples were taken out just before scanning and cleaned thoroughly to remove any traces of soil and dirt.The roots were scanned using an Epson Regent LA2400 scanner with a scale.The scans were analyzed using WinRHIZO (Regent Instruments Inc., Canada) for total root length (RL, mm.), diameter and root surface area (RA, mm 2 ).SRL (cm g −1 ) was calculated as the ratio of the total root length and the root's dry biomass.RMF (g g −1 ) was calculated as the fraction of root dry mass and total plant dry biomass.Seasonality index (SI) of each species was derived from a previous analysis that estimated species' preference to more of less seasonal conditions based on their probability of occurrence across a landscape-wide gradient of seasonal drought (Krishnadas et al., 2021).Positive value of seasonality index indicates that the species prefers less seasonal conditions, i.e., wet associated.

| Statistical analysis
We used mixed-effect models with appropriate error structures (see below) to test the impact of drought on survival, growth and photosynthesis.Preliminary models for growth showed that trends did not vary between growth rate ((final height -initial height) number of days between germination and harvest, Table A3) or relative growth rate (growth rate/initial height, Table A5) and we used the latter for interpretation since relative growth rates were standardized by initial height (See Table A4).
For the first question concerning differences in species response to drought in relation to their affinity for less vs. more seasonal areas, we modeled species' performance as a function of seasonality index, drought treatment, and their interaction.For question 2, we tested the role of each trait in governing drought response.We implemented separate models per trait, where species' performance was modeled as a function of their trait values, drought treatment, and drought-trait interaction.In all models, we used species-level trait means.Species ID and block ID were included as random intercepts.Survival was modeled using a Bernoulli error structure and Gaussian errors were assumed for growth and photosynthesis.
Because species differed in the extent of variation in their response, we added a weighting term to separately model the variance per species.Trait models for SRL were implemented with and without the species Actinodaphne malabarica, which had small seedlings with very light roots and high SRL values that were an outlier compared with other species.
Model structure in Wilkinson-Rogers notation: Finally, to understand the impact of the overall trait-based phenotype on performance during drought (question 3), we used a twostep approach.First, we performed a pairwise correlation among traits (Figure A2) and a principal component analysis (PCA) on traits to obtain composite phenotypes defined by trait combinations.
Then, we modeled individual survival, growth, and photosynthesis in control vs. drought plants using GLMM with species-level intercepts and slopes.From these, we extracted each species' baseline performance (survival, growth, photosynthesis) in well-watered conditions and the change in performance with drought and conducted a PCA on these six variables to get species' composite response to drought in relation to their baseline differences in growth, survival and photosynthesis.For both PCAs, we performed varimax rotation to simplify loadings for each factor.In the second step, to explore the correspondence between traits and performance, we performed Procrustean superimposition on the two PCA, using the axes that explained >75% of the variation in each PCA.Procrustes analysis checks the strength of association between two ordinations for a common object (species in our case), i.e., how well species positions in trait space matched their positions in performance space (Peres-neto & Jackson, 2001;Rüger et al., 2018).Significance of the Performance variable ∼ trait + treatment + trait: treatment + (1| species) + (1| block), weights = varIdent(form = ∼ 1| species) Procrustes correlation was assessed using permutation tests with 10,000 random iterations.

| Multi-trait phenotype and drought response
Ordination for species performance under drought, the first two PCA axes explained 77.2% of the variation in performance among species-level coefficients for survival, growth and photosynthesis (Figure A3 and Table A8).The first performance axis separated species according to their survival under drought vs. control conditions, with positive values indicating poor drought survival (Table A8).
As a whole, the PCA showed that species that grew well also had higher photosynthetic rates and had higher survival in both drought and well-watered conditions (Figure A3 and Table A8).
After varimax rotation, PCA for traits required four dimensions to explain ~80% of the variation among species (Figure 4a, Table A10).
Most axes had at least two traits with loadings >0.5.The first PCA axis separated species based on their stomatal, leaf structural and xylem traits and explained 29% variation.Positive values of axis 1 corresponded to lower SLA, larger stomates, higher leaf structural tissue, smaller xylem, and larger root mass fraction (Table S11).Axis 2 separated species based on water acquisition traits, explaining ~21% variation.Species with positive values on axis 2 had larger, longer and denser roots (Table A11).However, the third and fourth PCA axes reveal alternative combinations of measured traits.For instance, F I G U R E 3 Trait-mediated photosynthesis performance in drought.Seedling photosynthesis rates in drought vs. well-watered conditions as a function of traits.Photosynthesis rate was modeled using Gaussian error structure for seedlings under drought vs. wellwatered conditions as a function of (a) Xylem diameter, (b) Stomatal density, (c) Specific leaf area, (d) Leaf dry matter content, (e) Vein density, (f) Stomatal Size, (g) Stomatal length and (h) Root diameter with random intercepts for species and experimental block and weightage of species.Lines indicate predicted response and points denote raw data.Continuous lines indicate a significant relationship at alpha < .05 and dashed lines indicate non-significant results.Red lines and points denote drought conditions and blue denotes well-watered conditions.As. denotes the Photosynthesis rate.compared with axis 1, axis 3 comprised species with higher LDMC that also had higher SRL and smaller stomates.Axis 4 showed that, in contrast to axis 1, low SLA species could have smaller xylem and higher vein density (see Table A9).
Procrustes analysis between the ordinations of species performance (two axes) and traits (four axes) showed substantial concordance (sum of squares = 38.24, Figure 4b, r = .66,p = .03;permutation test for significance at alpha 0.05).This indicates that species' performance in drought conditions could be explained by composite phenotypes discerned from anatomical and morphological traits.

| DISCUSS ION
We found that drought reduced survival and growth for seedlings of 16 tree species from a wet tropical forest region in India (Western Ghats), but species varied substantially in their responses to drought.Interspecific variation in survival and growth with drought was primarily explained by anatomical traits that mediate water use and hydraulic safety.As expected, change in photosynthetic performance under drought stress correlated best with anatomical traits related to water use and gas exchange.Generally, species with fewer stomata and smaller xylem showed smaller declines in performance, growth and survival with drought.Growth constituted a distinct axis when compared with survival and photosynthesis for variation in species' drought response (Figure A3), but the overall variation in performance among species was explained by composite phenotypes determined by multiple traits.Species from more seasonal sites showed larger proportional decline in photosynthetic performance under drought (Figure 1c), but species affinity to seasonality conditions did not explain variation in drought response of growth or survival.
Different aspects of species' drought performance were regulated by different traits, but stomata and stem xylem traits correlated with all tested measures of drought performance, i.e., photosynthesis, growth and survival.Seedling survival and growth in well-watered conditions improved with larger xylem and stomata (Figure 1a,b).In well-watered conditions therefore, survival improves with traits that allow moving larger volumes of water for carbon assimilation.Accordingly, species with larger xylem also had higher photosynthetic rates (Figure 3a).The benefit of larger xylem for survival diminished in dry conditions (Figure 2a), consistent with larger xylem being more vulnerable to hydraulic failure due to embolism and cavitation, as known from mature trees (Hoffmann et al., 2011;Santiago et al., 2018;Skelton et al., 2015).
Tighter stomatal packing improved photosynthesis and growth when water was not limited (Medeiros et al., 2019).As with xylem size, this advantage reversed in drought conditions, when species with larger, more densely packed stomates experienced substantially larger declines in photosynthesis (Figure 3b,f,g) and growth (Figure 2e).Larger stomata imply higher conductance but also greater water loss, as with more numerous stomata (Lawson & Blatt, 2014), likely resulting in greater water loss, which becomes problematic during water limitation.Note that stomatal area fraction (SAF = SD*SS) showed results consistent with stomatal density (Figure A3), but smaller stomatal size did not correlate with higher stomatal density, suggesting that stomatal size and density form different axes to regulate water-use (Anderegg et al., 2018).
Species that proportionally invested in longer roots, i.e., lower specific root length, had better survival when water was plenty, but experienced larger decline in survival in drought conditions (Table A7i).This is consistent with resource-acquisitive being disadvantaged in stressful conditions, but this was not reflected in photosynthetic performance as expected.Photosynthetic rates did not vary with root mass fraction, or the length and area of roots in either well-watered or drought conditions, even as average photosynthetic rates declined with water scarcity.We thus found no support for the expectation that investing in below ground tissue can help maintain growth during drought by improving access to water (Chitra-Tarak et al., 2018;Joslin et al., 2000;Paz et al., 2015;Poorter & Markesteijn, 2008).It is possible that our findings may be an artifact of pot experiments where plants can only explore a finite soil volume, which might reduce the relevance of root traits in acquiring water.
In general, trait values that supported better growth in wellwatered conditions resulted in larger declines in growth with drought, suggesting a safety-efficiency trade-off with respect to water use (Choat et al., 2012).Species with larger xylem, suggesting higher hydraulic conductivity, also decreased photosynthesis during drought, but in this case, reduced photosynthesis translated to greater mortality suggesting carbon starvation with drought.
Variability of trait effects on different aspects of performance was unsurprising when considering that survival and photosynthesis comprised distinct axes when compared with the growth for variation among species (see PCA of species coefficients for performance variables, Figure A10).Moreover, different combinations of traits were possible and required to explain interspecific variation in the overall drought response represented by changes in survival, growth, and photosynthesis.
A few caveats temper interpretation of our results.Extrapolating from pot experiments to performance in field conditions merits caution.Plant size affects soil water depletion and transpiration rates when growing in a limited volume of soil, influencing comparisons among species (Comita & Engelbrecht, 2014).Initial height, however, did not predict growth rates or drought responses in our study.
Second, in closed canopy forests, shading affects air temperature and vapor pressure deficit (VPD), which can modulate the effects of water availability (Kupers et al., 2019).In addition, larger trees often transport water from deeper layers towards the surface (hydraulic redistribution) and increase water availability for seedlings to buffer drought stress.Alternatively, drought in natural settings may increase mortality due to other stressors such as nutrient limitation or herbivory, which we did not assess.Transplant experiments in the field will be needed to understand 'real world' implications of drought for seedling performance (Comita & Engelbrecht, 2014).

| CON CLUS ION
We found that anatomical traits related to xylem and stomata played a more prominent role than commonly used traits (e.g., SLA, LDMC, Greenwood et al., 2017) in explaining the drought response of tree seedlings from tropical humid forests of south Asia.
These anatomical traits are directly involved in water movement, gas exchange and hydraulic safety (Brodribb et al., 2003;Brodribb & Holbrook, 2003;Sterling, 2005), and may serve as effective proxies to link hydraulic performance and whole-plant response of tree seedlings to drought, allowing relatively rapid comparisons among species in field conditions and making generalizable predictions.Overall, drought conditions equalized performance differences seen in well-watered conditions or advantaged some species based on their traits.Increased frequency of intra-seasonal droughts or severity of the dry season due to global environmental change may thus influence which tree species regenerates successfully and accordingly alter functional composition of humid tropical forests (Chou et al., 2013).

ACK N OWLED G EM ENTS
Seeds from Sirsi were provided by Raghavendra H. V. and from Valparai were provided by Divya Mudappa, Srinivasan Kasinathan and others in the Nature Conservation Foundation team, for which we are grateful.We also thank the Greenhouse facility at The National Centre for Biological Sciences (NCBS) for the use of their infrastructure and personnel support.We thank Dr. Jayanti Ray Mukherjee from Azim Premji University, Bengaluru, for providing us the lab space and WinRHIZO for root trait measurements.We thank Abhishek Gopal for help with initial data analysis.Funding for this project came from the NCBS campus fellowship (through TIFR and the Dept. of Atomic Energy, Government of India) and National Geographic Society explorer level-I grant.We acknowledge the support provided by NCBS, TIFR and the Dept. of Atomic Energy, Government of India.

FU N D I N G I N FO R M ATI O N
NCBS campus fellowship (through TIFR and the Dept. of Atomic Energy, Government of India); National Geographic Society explorer level-I grant.

CO N FLI C T O F I NTER E S T S TATEM ENT
The authors declare no conflict of interest.

O PE N R E S E A RCH BA D G E S
This article has earned an Open Data badge for making publicly available the digitally-shareable data necessary to reproduce the reported results.The data is available at doi: 10.5061/ dryad.4j0zpc8n1.

S U PP O RTI N G I N FO R M ATI O N
Additional supporting information can be found online in the Supporting Information section at the end of this article.
How to cite this article: Jhaveri, R., Cannanbilla, L., Bhat, K. S. Note: Species code is the combination of first three letters of the Genus followed by first three letters of the Species.These codes are used to refer to each species in the main text.Higher seasonality index values indicate species associated to more seasonal forest (Krishnadas et al., 2021).Seeds used for the experiment were collected from two locations, Sirsi (Karnataka state, India) and Valparai (Tamil Nadu state, India).

TA B L E A 2
List of traits quantified in the study and their abbreviations as used in the manuscript.SI values were taken from Krishnadas et al. (2021).Intensity of the color corresponds to the strength of r.Abbreviations for the traits are given in Table A2.

F
Trait-mediated growth rate and survival in drought.Seedling survival in drought vs. well-watered conditions as a function of (a) Xylem diameter, (b) Stomatal size and (c) Leaf dry matter content.Survival was modeled using Bernoulli errors with random intercepts for species and experimental blocks.Seedling growth under drought vs. wellwatered conditions as a function of (d) Xylem diameter, (e) Stomatal density, (f) Leaf dry matter content, (g) Vein density and (h) Root diameter.Growth rate was modeled using Gaussian errors with random intercepts for species and experimental block and weightage of species.Lines indicate predicted response and points denote raw data.Continuous lines indicate significant relationship at alpha <.05 and dashed lines indicate nonsignificant results.Red lines and points denote drought and blue denotes wellwatered conditions.

F
Comparing stomatal trait-mediated seedling performance in drought.Seedling performance in drought vs. well-watered conditions as a function of stomatal density, stomatal size and stomatal area fraction.Seedling survival under drought vs wel watered condition was modeled using Bernoulli errors as a function of (a) stomatal density, (b) stomatal size and (c) stomatal area fraction.Growth rate was modeled using Gaussian error structure for seedlings under drought vs. well-watered conditions as a function of (d) stomatal density, (e) stomatal size and (f) stomatal area fraction.Photosynthesis rate was modeled using Gaussian error structure for seedlings under drought vs. well-watered conditions as a function of (g) stomatal density, (h) stomatal size and (i) stomatal area fraction.For all the models, random intercepts for species and experimental block were used with and weightage of species.Lines indicate predicted response and points denote raw data.Continuous lines indicate a significant relationship at alpha <.05 and dashed lines indicate non-significant results.Red lines and points denote drought conditions and blue denotes well-watered conditions.F I G U R E A 2 Plant trait correlation matrix: Matrix showing correlation coefficient (r) among plant traits.Color scale corresponding to the strength and sign of r is given on the right.Red denotes negative and blue denotes positive r.

F
Seedling performance and survival in drought and well watered conditions.Principal Component Analysis decomposing variation among seedling performance and survival in drought and well watered condition.Direction of the arrows indicate the relative influence and the length of the arrow indicates the magnitude of influence of each predictor on seedling performance and survival.Survival indicates modeled survival probability, Growth indicates the modeled Relative Growth Rate and As.indicates modeled Photosynthetic rate.Suffix for each predictor indicates well watered (C) or drought (D) treatment.

Growth Rate ~ SD * Treatment (Drought/Control) Conditional R Squared = .2089 β Std. Error t-Value p-Value DF
Trait-mediated growth rate in drought.

Growth Rate ~ RF * Treatment (Drought/Control) Conditional R Squared = .3712 β Std. Error t-Value p-Value DF
List of model coefficients for seedling growth rate under drought vs. well-watered conditions as a function of traits.Sub-tables a-l refer to different models with each of 12 traits (refer to TableA2for trait abbreviation).Growth rate was modeled using Gaussian errors with random intercepts for species and experimental block and weightage of species.First row of each sub-table lists model response and explanatory variables with appropriate relationship.Model Conditional R-squares are given in the second row of each sub-table.Model co-efficients (β), standard error (Std.Error), t-statistics (t-value), p-value and degree of freedom (DF) are listed in columns in rows corresponding to each explanatory variable.Seasonal affiliation-mediated performance and survival in drought.Seedling survival and performance in drought vs well-watered conditions are modeled as a function of their seasonal affiliation.(a) Results from generalized linear mixed effects model with binomial error for seedling survival in drought vs well-watered conditions as a function of their seasonal affiliations with random intercepts for species and experimental block.Results from linear mixed effects model for seedling relative growth rates (b) and photosynthesis rate (c) in drought vs well-watered conditions as a function of their seasonal affiliation with random intercepts for species and experimental block and weightage of their seasonal index.First row of each sub-table lists model response and explanatory variables with appropriate relationship.Model Conditional R-squares are given in the second row of each sub-table.Model co-efficients (β), standard error (Std.Trait-mediated relative growth rate in drought. Error), t-statistics (t-value), z-statistics (z Value), p-value and degree of freedom (DF) are listed in columns in rows corresponding to each explanatory variable.TA B L E A 5

Relative Growth Rate ~ SD * Treatment (Drought/Control) Conditional R Squared = .28 β Std. Error t-Value p-Value DF
Trait-mediated Seedling Survival in drought.

VD * Treatment (Drought/Control) Conditional R Squared = .32 β Std. Error z Value p Value
Two Principal Component Analysis dimensions with three performance parameters in drought and well watered condition, explained ~77% variation among species.Suffix for each predictor indicates well watered (C) or drought (D) treatment.Species vice PCA loadings for performance model and survival model coefficients.Two Principal Component Analysis dimensions with 16 species, explained ~77% variation among species performance in drought vs. well watered conditions.PCA loadings for species traits.Four Principal Component Analysis dimensions with 12 traits explaining ~80% variation among species traits.Species vice PCA loadings for traits.Four Principal Component Analysis dimensions with 16 species, explained ~80% variation among species traits.
TA B L E A 1 0 TA B L E A 11